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SYIVOPSIS 


This work considers the prohlem of detormining the 
maximal arc disjoint flow in a multicommoditj network with 
one source node for each commodity and a common terminal 
node for all the commodities. This prohlem is a special case 
of the general multicommodity flow prohlem and considers the 
restriction that more than one commodity should not flow on any 
arc. The problem is essentially that of partitioning the 
arcs of the network into as many sets as the number of commo- 
dities and allowing only one commodity to flow over each sot, 
such that the sum of flows of all commodities from their 
source nodes to the common terminal is maximum. 

The mathematical formulation of the prohlem is done as a 
mixed integer linear program. Two algorithms based on the well 
known branch and hound strategy are proposed for the solution 
of the problem. The two algorithms differ from each other only 
in the method of determining the upper hound of each node in 
the branch and bound tree. 

Computer codes are written for the proposed algorithms in 
PORTIhm-lO and implemented on DEG-1090 time sharing computer 
system. Problems of size as large as 200 nodes, 300 arcs and 
5 commodities were solved. Even fox the largest size— problems 



solved, the avor.-^ge CPU time -was less than 6 seconds in 
case of both the algorithms^ The computational performances 
of the proposed algorithms for a large number of randomly 
generated net'works of varied sizes are presented. 
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II^RODUCTIOF 

An important area in the study of netw6<rks is -that of the 
maximum possible flow between two specified nodes of the network. 
In the broad sense, a flow on a graph is a way of sending 
objects from one vertex to another by travelling along the arcs 
in their directions. 'The problem of finding the maximum flow 
can be posed in almost every instance in which commodities flow 
from a source node to a destination node. Thus in traffic net- 
work, we may be interested in the maximum rate of traffic flow 
between two cities, in an electric power distribution netv^ork, 

■we may be interested in the maximum power transmitted from a 
generating station or stations to' a particular location; in a 
natural gas distribution netv/ork we may be interested in the 
maximum rate at v/hich gas can be supplied to a particular 
Consumer etc. 

This is the most basic problem in network flov; ^'fhose 
solution is readily available through a labelling procedure 
due to lord and lulkerson [51 or through the method of 
preflows by Karzanov 171. There are two important extensions 
to the one source - one terminal flow problem. 

The first is the multiterminal flow problem in which 
several nodes are designated as source-terminal (s,t) pairs. 
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vith the same commodity flovs through the netx«?ork. Por example, 

I 

this is the case of a telephone network^ in which any unordeped 
pair of n cities covered by the net\A.'ork may indeed serve as 
the (s-t) pair. 

The second is the 'multicommodity' flow problems in v;hich 
several commodities flov? simultaneously from designated sources 
to designated terminals. The sources and the terminals may be 
different for different commodities. The problem is to find such 
a flow that the sum of the different flows from the respective 
sources to the respective terminals is maximum. For example, 
in multiproduct production - distribution systems, the produ- 
ction units and the customers serve as the source-terminal pairs 
and the distribution channels, the arcs of the network. 

The problem considered in this work is a variation of the 
multicommodity flow problem where all the commodities flow to 
a common terminal node. Here an additional constraint that the 
flows are to be arc-disjoint is also imposed. This means no 
arc allows the flow of more than one commodity over it or each 
commodity will flow over mutually exclusive set of arcs. For 
example, in pipe line networks, where more than one commodity 
flow and these commodities originate from different sources 
and terminate in a common processing unit, there could be arc 
disjoint constraints due to technological reasons. 
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The sarae problem can also be posed as the problem of 
partioning 'the set of arcs into as many mutually exclusive 
sets as the number of commodities that are to pass through 
the network. Each of these subsets is assigned to one commo- 
dity. On such a subset, the maximal flow of the commoditj^ for 
vjhich it is assigned could be found using any sino-le commodity, 
single source, single terminal maximal flow algorithm. Thus, 
it is required to find that partition which gives the maximum 
sum of flows. 

¥e could also think of two other similar maximal flow 
problems in a network. One is the maximal arc disjoint flow 
problem in a multicommodity common source multiterminal net- 
work. The other is the maximal node disjoint flovj' problem in a 
multicommodity common terminal (or common source) network in 
which more than one commodity is not allowed to pass over any 
node (other than the common tejcminal (common source) node) in the 
network. Because of the underlying similarities, these problems 
could also be solved using the solution procedures suggested for 
maximal arc-disjoint flow problem in a multicommodity, common 
teminal network. 

The common source problem could be solved using a branch 
and bound method similar to the one suggested for common ter- 
minal problem. Instead of connecting all the source nodes to a 
supersource, in the case of common terminal flow problem, all 
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the different sink nodes are connected to fictitious super- 
sink using arcs of infinite capacity. 

The prohlem of node-disjoint flov.' is a further restriction 
of the arc disjoint flov prohlein. That is, all node disjoint 
flows are arc disjoint too. (But the converse of this statement 
that all arc disjoint flows are node disjoint is not true). 
Moreover all node disjoint flow problems may be converted to 
arc disjoint flow problems by a simple transformation given 
below. 

For each node K (other than the source nodes and the 
common terminal node), add one new node S'. K and S' are 
connected by a directed arc of large capacity. All the inco- 
ming arcs to the corresponding node in the original network 
are incident upon the node K and all the outgoing arcs from 
the corresponding node in the original netv/ork will emerge out 
from the node S' . The source nodes and the common terminal 
node are kept undisturbed. Thus, if the original network had 
n nodes, m arcs and p commodities, the transformed network will 
have 2n-p-1 nodes and m+n-p-1 arcs. This transformation is 
illustrated in Fig. 1. 

On the transformed network, the maximal arc disjoint flow 
will be the same as the node disjoint flow in the original 
network. This becomes obvious if we consider the flow over the 
newly added arcs in the transformed network to be the flow 
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Newly added arcs of infinite 

capacity 


Pig. 1 I'ransf ormat ion of node-dis;ioint flow problem to arc- 
disjoint flow problem 

through the node in the original network corresponding to that 
arc in the transformed network. Since the flows -oltained are arc 
disjoint there is only one commodity flowing over the newly 
added arcs in the transformed network w^hich makes the flow to he 
node disjoint in the original network. Since the newly added 
arcs are given infinite (a large number) capacity, nodes do not 
restrict the amount of flow. But solving the problem on the 
transformed network suffers from the drawback that the transformed 
network has more arcs and nodes than the original network. The 
problem of arc disjoint flow could also be solved using a 
branch and bound method similar to the one suggested in this 
v;ork for arc disjoint flows. At each node of the branch and 
bound tree, instead of assigning an arc to a commodity, a node 
is assigned to a commodity. Obviously this will provide a 
solution to the problem with lesser amount of computation than 
the one provided by applying the solution procedure for arc 
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disjoint flow in the modified network, because the number of 
arcs and nodes in the original network is lesser than the 
number of arcs and nodes in the modified network. 

The solution procedures for maximal arc disjoint flow 
problem in multicommodity common source network and maximal node 
disjoint flow problem in multicommodity common terminal (or 
Common source) network are not discussed any further in this 
work as they are similar to the solution procedures for the 
maximal arc disjoint flow problem in multicommodity common 
terminal network. 
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IITERilTURB SURVEY 

The class of netv^ork optimization problems, called the 
multicommodity flow problems arise in network modelling wherever 
commodities, vehicles or messages are to be shipped or transmi- 
tted from certain source nodes to certain terminal nodes of the 
underlying network. Recent applications of mathematical pro- 
gramming techniques to traffic equilibrium problems in trans- 
portation studies as well as computer networks analysis has 
generated considerable interest in this class of problems. 

The multicommodity maximal flow problem arose out of the 
single commodity)" version. The efficiency of the labelling 
method for solving the single commodity maximal flow problem, and 
the theoretical insight offered by the Max-Elow-Min-Out theorem 
led investigators to seek natural generalizations to the multi- 
commodity version. Unfortunately such generalizations do not 
preserve the nice combinatorial features of the single commodity 
problem. 

A multicommodity cut or disconnecting set is a set of 
arcs upon the removal of which no chains exist from the source 
to the terminal for any commodity. The multicommodity min-cut 
is one with minimal capacity. Because of weak duality, the 
maximal flow is always less than or equal to Minimal-cut and 
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the netvork is called '’gapless" if equality holds. Single 
Commodity problems are gapless by the Max-flo-w-Hin-Cut theorem, 
but the inequality can be strict for multicommodity flow 
problems. The multicommodity flow problem is not necessarily 
integer either. This means that the multicommodity flow 
problem may fail to have integral optimal solutions although 
arc capacities are integral. The integrality property for the 
single commodity problem followed from the total modularity 
of the constraint matrix. Jewell [6] has tried to extend this 
property to mult i commodity flows but has met little success. 

However, exploiting the combinatorial structure of the 
multicommoditj?' problems, a few special results have been 
developed. Host of the special results are restricted to 
special networks or networks with ntmost two sources or two or 
three commodities. Using matroid arguments, Evans et al. ( 2] 
show that a multicommodity transportation problem has a uni- 
modular constraint matrix if and only if the number of sources 
or sinks is not greater than two. 

The two commodity max flow problems has been investigated 
extensively^ Algorithms for finding the maximal two-commodity 
flow on undirected graphs arc given in Hu [51 . Hu shows that 
this problems has integral solution for even arc capacities. It 
is hard to extend these results to problems with three or more 
commodities except in very special cases, for three-commodity 
flow problem, Rothfarb et al [12J prove that a special three- 
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coiTimodity undirected graph in v/hich all nodes are either 
sources or sinks is gapless and give an algorithm, to find the 
maxinial flow. 

Algorithms are available for other specialised maximal 
flow prohleuis. Kleitman [S] solves the case where each node 
is a sink for all commoditios hut one. Rothfarah et al [11] 
consider the multicommodity network flow rrohlem where oil the 
Commodities are destined for a common terminal. Thev have 
presented an, algorithm to maximize J ? for decreasing 

2 k 

weights a > a > ... >« ^ where K is the number of 

commodities. 

Ford and Pulkerson [41 have suggested solution techniques 
to the multicoamiodity maximal flow problems. They proposed 
the arc-chain formulation of the problem, and solved the 
problems using revised simplex method. 

' The problem of arc disjoint flow in multicommodity 
network was studied by Ahuja [1 I . He suggested a branch and 
bound algorithm for the solution of this problem. The upper 
bound was obtained by solving p independent single commodity 
maximal flow problems where p is the number of commodities. 

A thorough survey of the literature suggests that 
practically no work has been published in the area of arc- 
disjoint flow in multicommodity networks with a common terminal 
node. This acted as the prime motivation for the present work. 
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M. JT HIS'EiT I OiiE EOffi-IUM lOR 

Most of the network flov; prohleras lend themselves in a 
netural manner to linear Programming formulation. Therefore, 
it is possible to solve these problems using the available 
techniques for the solution of linear programming oroblems. 
Hov’ever, due to the special structure of the network flow 
problems number of efficient algorithms and elegant theorems 
have been developed. In most of the network flow problems, 
the optimum solutions are always integers, a fact which is not 
true in linear programs. But this property of the single 
commodity flow that the optimal solutions are integers is not 
always true for multicommodity networks. Yet it is to be 
noted that the multicommodity flow problem taken up in this 
work has integral solutions. This is clear from the fact that 
the solution to the problem is the sum of the solution of 
p single commodity single source, single terminal flow 
problems (p is the number of commodities to be flown in the 
network) solved over p networks, each of which is made up 
of mutually exclusive set of arcs. Solutions for each of these 
individual problems are integral and so the sum of p inte- 
gral solution will also be integral. Exploiting this special 


property of the linear program solution, we could develop better 
solution procedures than the linear program formulation of the 
problem. 
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The formulation of the problem as a mixed integer linear 
program is given belovj. 

The network: GCi?,!) is a finite set of nodes I'J’ containing 

n elements and a finite set of ordered pair of nodes called 

arcs, A containing m ‘elements, l-'fe shall use to indicate the 

ith node and A. . to indicate the directed arc leading from 

J 

U. to IT.. A network is connected if for every partitioning of 

X J 

nodes, of the network into subsets Z and X, there is either an 
arc A. . or A., with IT. eX and iT. e X. The word 'network' will 

-L J J -A J 

be used to denote connected network throughout this v/ork. It 
is also assumed that all the arcs are directed and. there is 
atmost one directed arc between any pair of nodes. That is 
both A. . and A., will not be members of the set of arcs. 

1 J J X 

Clearly, the problem of maximum flow in a network is meaning- 
ful only if the arcs of the .network or a subset of it possess 
upper limits called the capacities of the arc. let each arc 
A. . have a capacity b. . which is the maximum amount that can 

X J 1 J , 

flow over that arc. 

In the multicommodity network considered, more than one 
distinct commodity flov/ over the network with as the origi- 
nating node (source node) for the ith commodity and 't' as the 
destination node (terminal node) for all the commodities. Let 
p be the number of commodities that are to flov; over the net- 
work. So there will be p source nodes and one common terminal 
node in the network. 
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Associated vith each aode of the netvork, there are two 
integer numbers called in-degree and out-degrce. In -degree 
of a node is the total number of arcs coining into that node 
and the out-degree is the total number of arCvS going out of 
that node. It is also assumed that the in-degree of all the 
source nodes and the out-degree of the common terminal node are 
zero. 


The problem can be posed as a mixed integer linear 
programming problem. 

Let denote the flow of the kth commodity between the 
source node for that commodity and the common terminal. Let 
x^^ denote the amount of flow of the kth commodity over the arc 




— Ic 

Let us introduce the zero-one variables y. . 

J 


wh ere 


y^^ = 1 if the kth commodity passes over arc A. . 

X J X J 

= 0 otherv^ise. 


' k 

There will be mxp variables of the type x. - and another 

X J 

k 

mxp variables of the type y- • where m and p are the number of 

X J 

arcs and the number of commodities respectively. 


Ob,i e ctive f unc tion : 

The objective is to maximize the sura of the flows of all 
commodities. So the objective function will be 

^ k 

Maximize v 

k=1 
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Oonstraint s 

The constraints of this prohlen are of 4 types. They are 
given helow. 

1. Plow conservation constraints : 

These constraints represent the fact that for each 
commodity, the total flow emerging out of its source node and 
the total flow terminating in the common terminal node are 
numerically equal to the total flow of that commodity and for 
all nodes other than the source nodes and the common terminal 
node, the flov? is conserved (i.e. the total incoming flow is 
equal to the total out-going flow). 


These constraints can be expressed mathematically as 


n n 

f „.k 
f V 

if 

i 

= s. 



k 

I xj, - I 

X . . = k 

n i-v 

if 

i 

= t 

3=1 3=1 

i 

1 0 

if 

i 

is any node other than 


the source nodes and the common 
terminal node, 
for k = 1,2, ,,, p, 

Corresponding to each node of the network, there are p 
constraints of this type. Thus the total number of flow 
conservation constraints is nxp. 
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2. Ixro disgoiiit constraiats ; 

These constraints represent the restriction that more 
than one commoditj should not flov/ over an arc. This restri- 
ction can be expressed mathematically as 

P 

I f 1 for all(i,3) e A, 

k=1 

On the arc 1-^, flow of the kth commodity can be positive 

X J 

k 

only if the value of the variable y - • is 1. for any other 

X J 

commodity k’ (i.e. k' = 1,2, ... p, k) because of the above 

k ' 

constraint, the corresponding zero-one variables y. . takes 

X J 

zero value; thus indicating that there is no flow of the 
commodity k' over the arc 

X J 

for each arc in the network, there is a corresponding 
arc disjoint constraint. Thus the total number of arc disjoint 
constraints is m. 

5. Capacity constraints : 

These constraints ensure that the total flow over an arc 
will not be more than the capacity of that arc. 

Since the arc disjoint constraints take care that only 
one commodity flows over an arc, the capacity constraints 
have to only restrict the flow of any commodity over the arc 
A., to be not greater than its capacity b. (It is not nece- 

X J X J 

ssary to restrict the sum of the flows of different commodities 
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over the ^.rc to he not greater than the capacity of the arc). 
The matlieniatioal expression of these constraints is 

k; = 1,2, . . . , p . 

Corresponding to each arc in the netv^orlc, there vdll he 
p such constraints. Therefore, the total nuaiher of cappcity 
constraints are mxp. 

4. ilon-negat ivity constraints ; 

The non-negativity constraints restrict the flow of any 
commodity over an arc to he non-negat ive. Thus, 

> 0 for all (i,g) e A 

k = 1,2, . , . , p 

fhere are mxp constraints of this t^^pe. 

Also 

= 0 or 1. 

There are another mxp constraints of this type too. 


ns seen in the formulation, this is a mixed integer linear 

Ic 

programming problem where the variables y. • are restricted to 

1 J 

0 or 1. ^or a network v/ith m arcs, n nodes and p commodities, 


the total number of variables in the mathematical formulation 


will be 2xmxp. Out of this mxp are zero -one variables. 
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There will be (Janxp) + (n:q)) + ra constraints of which mxp 

Ic 

constraints are the constraints v/hich restrict x. ■ to be non- 

i. J 

negative and another mxp constraints are the constraints which 
Ic 

restrict to take values either 0 or 1 . 

X J 


The problem can be solved by any one of the mixed integer 
programming algorithms. However, the solution of a moderately 
sized problem by such an ax^proach may be prohibitive because of 
large computer raemory requirements and computational complexity. 
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SCIUTIOIT METHODOLOGY 

In this •'hapt er, two algorithms are developed for 
solving the maximal arc disjoint flow problem in a multi- 
commodity Common terminal network. Both the algorithms are 
based on the well known branch and bound strategy. The 
branch jmd bound approach is a veiy useful tool for solving 
combinatorial optimisation problems. In the process of 
seeking the optimal solxition, it iranlicitly entimerates all 
possible solutions. It performs it systematically and aims 
to do it with least computations. 

4.1 Branch and Bound Philosophy 

In the course of applying the branch and bound approach, 
the overall set of feasible solutions to the problem, is 
partitioned into many simpler subsets. At each stage, one 
promising subset is chosen and an effort is made to find the 
best feasible solution from it. If it is found, then that 
subset is said to be fathomed. If the best feasible solution 
is not found, then that subset is again par'citioned into two 
or more simpler subsets (this operation is knov/n as branching) 
and the same process is repeated again. At each partitioning 
of a set into its subsets, a lower bound (or upper bound) 
which gives the minimum value (or maximum value) of the 
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ob.iective function vjhich any feasible solution in that 
subset may yield. If this bound is greater than (or 

lov;er than) some best known feasible solution (known as 
incumbent) then that subset is not further branched. This 
operation is known as pruning. When all the subsets are 
either fathomed or pruned, incumbent provides the value of 
optimal solution. 

4.2 Algorithm - 1 

4.2.1 Upper bounding .strategy 

Let us consider a netv;ork through which p commodities 
are to flow. At any node of the branch and bound tree 
developed for the network, the arcs of the network can be 
dividied into (p+1 ) groups. Let A^ be the set of arcs which 
are assigned to no commodity, further let A^ , A 2 , A^ , ... 

Ap represent the set of arcs assigned to commodities 1,2,3, ... 
p, respectively. 

It is to be noted that the sets A^ , A 2 , ... A , A are 

ir 

mutually exclusive and collectively exhaustive sets. 

To find the upper bound for each node of the branch 
and bound tree, a super source is added to the network and 
each of the source nodes are connected to this super source 
with arcs of infinite (large) capacity. The maximal flow 
from the super source to the common terminal which could he 
found out hy a modified maximal flow labelling algorithm in 
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which commodity 'i 'is allov/ed to pass only through the set 
of arcs is used as the upper hound flow for the 

suhprohlem corresponding to that node of the hranch and 
hound tree. 

It is to he noted that for finding the upper hound of 
any suhprohlem, we are not solving p, single commodity 
single source, single terminal flows independently v^hich 
would involve the following procedure. 

I'or commodity 'i', remove all the arcs which are 
assigned to other commodities from the netv;orte and find the 
maximal flow from the source of the ith commodity to common- 
terminal. The same process is repeated for all values of 
'i'. The sum of the p flow values thus obtained will logi- 
cally form an upper hound for the problem. But in this case, 
both the capacity constraints and arc disjoint constraints 
arc violated and so v/ill provide only a loose upper hound 
for the problem. 

The upper hound obtained by the proposed strategy will 
he close to the solution, as in the flow corresponding to 
this upper hound, only the arc disjoint constraints are 
violated for the arcs in the set hut the capacity con- 
straints are maintained for all arcs. This tightness of 
the upper hound accounts for the efficiency of the branch 
and hound algorithm suggested. 
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Logically the upper hound value decreases in further 
assignment s , as at each stage an unassigned arc get 
assigned to one of tlie commodities. 3o its flexibility of 
getting assigned to any of the commodities is lost. 

The proof for the validity of this upper bound for 
the maximal flow for any particular assignment (i.e. at any 
particular node of the branch and bound tree) is given below. 

Proof ; 

Obviously the labelling algorithm terminates when no 
augmenting path exists between the source of ith commodity 
and the common terminal over the set of arcs for all 

values of 'i'. Let v be the total flow value of such a flow 
obtained. 

Let there exists another flow for the same assignment 
of arcs for which the total flow value is v^ and v^ > v. 

This means that for at least one of the commodities, 
the flow from its source to the common terminal should bo 
more than the present flow. But, since the algorithm ter- 
minates only v/hen no augmenting path is possible from any 
source to the terminal, an increased flow can be achieved 
only if the flow of some other commodities are reduced over 
arcs vdth which these commodities share at least one arc. 

Very obviously, these arcs Cronnot be in the sets ... A^ 
as these arcs allow the flow of only the commodity to which it 
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is assigned. So to get an increased flow of some commodities 

from their source nodes to the common terminal node, flow of 

those commodities over the arcs in set A are to he increased 

o 

and flow of other commodities over the arcs in set are to 
he decreased and the increment in flow obtained should he 
greater than the decrement of flow of other commodities. 

Since all commodities are allowed to flov; over the 
arcs in set A^, if any such augmenting paths had been 
possible, the labelling algorithm v/ould have detected it 
before termination, -it termination of the labelling algori- 
thm no such augmentations are possible and so there con not 
exist a flow for which the total flow value, v-i is greater 
than V, for the sarae assignment of arcs to different commo- 
dities, 

4.2.2 Branching Strategy 

Branching is done at each node of the branch and bound 
tree, considering the flow over all the arcs in the solution 
of the maximal flow problem at that node. Of all the arcs 
which allows simultaneous flo^'f of more than one commodity, 
the arc which has the maximum total flow (i.e., the sum of 
flows of all commodities) is choosen for branching). The 
arc thus selected serves as the basis of branching at that 
node and as many new nodes are created in the branch and 
bound tree as the number of commodities, for each node 
thus created, the arc selected for branching is assigned 
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seperatelj to cacli of the commodities. 

After e'^ch br^'nohing, the newly generated nodes of the 
branch and bound tree are examined and the feasible solu- 
tions are used for updating the incumbent. 

4.2,3 Searching Strategy 

The searching strategy selects the node of the branch 
and bound tree at which branching is done. Of the nodes 
■which are added latest to branch and bound tree, the one 
which has the largest upper bound ip chooser for further 
branching. This results in what is commonly termed as 
'depth first search strategy'. In this method, invstead of 
scanning every node which is directly connected to (adjacent 
to) the current node, we move to one adjacent node, as soon 
as possible , leaving the current node with possibly unexplored 
adjacent nodes for the time being, in other words, we trace 
a v^alk through the branch and bound tree, going on to a new 
node whenever possible. 


4.2,4 Notations Used : 


A 


P 

S 

S 


T 

^ij 


set of arcs in the network 

- number of commodities to be flown in the network 
source of the ith commodity 

- super source added to the network 

- common terminal node in the network 

- flow of the kth commodity over arc (i-j) 
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L 


I 

U(Pi) 


The set of arcs assigned to the ith comnioditj 
at any node of the branch and bound tree 
the set of arcs v>’hich are -^issignecl to none of 
the coamodities at any node of the branch and 
bound tree 

- }- set of active infeasible sub problems (p^) 
>7hich form the unfathomed nodes of tlio branch and 
hound tree 

- value of the current incumbent 

- upper bound of the subproblen p. 

X 


4.2.5 Step-wise Description of the iilgorithin 


Step - 0 ! 


Add a super source B and connect each source node to 
this super source with arcs of infinite capacity. 

Set 1=0 

Set I = {0} 

Set .1^ = {All arcs in the network) 

Set A^ ={0} for i = 1,2, ... p. 

Set = 0 for all (i,o)eA, 

k = 1,2, ... p. 

Solve the maximal flow problem from super source S 
to common terminal T. If the flow values are feasible, i.e. 
no arc has more than one commodity flowing over it; STOP 
(This is the optimum flow). Otherwise, let represent this 
subproblem, Md Pq 1:0 the list L. 
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Step - 1 : 


Check LiF?t L. If L = {0} , go ‘^o Step 6. Otherv;ise, 
tetke out the problen which has the largest upper hound of 
all the problems which are added latest to the branch and 
bound tree. Let it be p^^^. Update list L. 

Step 2 • 


If U(P^) < then prune that node of the branch ^nd 

bound tree and go to Step - 1. Otherwise scan all (i,;i)e4. 


Determine tho arc which flov; more than one coaciodity over it 

^ k 

^ •among all such arcs. 

k=1 


and has maxinum 


Step - 3 : 

Generate new vsubprobleras (new nodes for broach and bound 
tree), p^^^ , p^p » ••• Pnip* Gubproblem p^^^^ has the arc selected 
for brranching assigned to commodity 'i’. 

Step - 4 : 

Solve the subproblems Pin2> *•* Pmp subproblems 

corresponding to the newly generated nodes of the branch and 
bound tree) for maximal flow. Commodity ’i’ is allowed to 
flow only over the sot of arcs u The maximal flow 

thus obtained for the subproblem p . is used ag U(p .)• 

IXX J** 

Check whether the problem has feasible flow. If so; go to 


Step 5 
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j.xdd. all the sub p roll ecns vdth infeasible solutions 
to the list 1 in the increasing order of their corresponding 
upper bounds U(p .). Go to Stop - 1. 

step - 5 : 

If U(p .) >1, set I = U(p .)• G-o to Step - 4. 

‘-iTLi ’ ni 

Step - 6 : 

Current incumbent is the maxiiaal flov; value. Stop. 
4.2.6 Example 

I'he suggested procedure is illustrated with the help of 
0 small network; shov;n in Eig. 2. The three commodity network 
considered has 10 nodes and 12 arcs.411 the source nodes are 
Connected to a super source v/ith fictitious arcs of infinite 
(a large number) capacity shown by dotted lines. The cir- 
cled numbers placed near the end of the arc is the arc 
number and. the uncircled number placed near the beginning 
of the arc is the capacity of that arc. 

The v/orking of the nlgoritho is described step-wise. 

Step - 0 : 

Set 1=0 and I to tbe empty. All arcs are included 
in the set A^. and A^ are set to be empty. . The 

maximal flov; from the super source to common terminal found 
out using labelling algorithm is 25. Plow values are 
infeas,ible, let the current assignment of arcs correvspond 
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to subproblera (1) vuth upoer tound v^lue of 25. This 


prohlem is added to List L. 



Lig. 2 An illustrative example 


Stop - 1 ; 

The List L is checlced. It contains one subprohlom and 
it is taken out. List L is updated. 

Step - 2 ; 

Since upper hound of the suhprohleni is greater than 
the incumhont, all arcs are scanned and arc 11 is choosen 
as the basis of branching because it is the arc which allows 
the maximum total flow of all the arcs which allow the flow 
of more than one commodity. 

Step - 3 : 

The nodes 2,3 and 4 are generated in the branch and 
bound tree. Sodes 2,5 and 4 correspond to subproblems in 
which arc 11 is assigned to commodities 1,2 and 3 respectively. 
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Step -- 4 : 

Solving the suhprobloQS corresponding to nodes 2,3 nnd 
4, it is found thnt U(2) = 24, U(3) = 16, U(4) = 23. Sub- 
problem corr^L-sponding to node 3 has feasible solution. 
Therefore, go to Step - 5. 

Step - 5 : 

Since U(2) > I, I is set to U(2) = 16. Go to Step - 4. 

Step - 4 ; 

The infeasible subprobleras are added to LivSt L in the 
increasing order of their upper bounds. Thus subproblems 
corresponding node 2 becomes the I'^'st element in List L. Go 
to Step 1 . 

Iteration - 2 ; 

Step - 1 : 

The last element of list 1 (i.e. the subproblem 
corresponding to node 2) is taken out from List 1, 

Step - 2 : 

Since the upper bound of this subproblem is greater 
than the incumbent, scan all the arcs and arc number 12 is 
choosen to be the basis for further branching. 

Step - 3 

New nodes 5,6 and 7 are created in the branch and bound 

tree. 
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Step - 4 : 

The Golutioa ox the subprohlef/is j^'lelded the xollovdng 
values for upper hound U(5) = 15, U(6) =19, U(7) = 21. 

Since all the problems have feasible flov;, go to Step - 5. 

Step - 5 : 

Since U(5) < I, node 5 is pruned, since U(6) > T, 

I is set to U(6) = 19. Since lJ( 7) > I, I is set to U(7) = 2 
Go to Step - 4. 

Step - 4 ; 

Since none of the subproblens corresponding to the 
newly generated branch and bound nodes are infeasible, 
nothing is added to the List L. 

Iteration - 3 ; 

Step - 1 : 

Last olement in livSt L (i.e. the subproblem corres- 
ponding to node 4 of the branch and bound tree) is taken 
out . 

Step - 2 : 

Upper bound of this subproblem is greater than 
incumbent, .irc 12 is Selected as a basis for branching. 

Step - 3 : 

3 new nodes 8,9 and 10 are created in the branch and 


boimd tree. 
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Step " ^ : 

The p-uh^^rotlons corresponding to the nev^^l3^ generated 
branch and bound nodes are solved. U(8) = 16, U(9) = 20, 

TJ(10) = 19. Since all these subprobleras have feasible- 
solutions, go to Step “ 5. 

Step - 5 ; 

The upper bound for all the three subproblens are less 
than the current incuvPoent and -^o all those nodes vore 
pruned. Go to Step - 4. 

St ep - 4 5 

ilone of the newlj generated subproblems gave infeasible 
solutions. So nothing is added to list I. Go to Step 1. 

Iteration - 4 : 

Checked the list I and it is found empty. Stop. 

The branch ana bound tree developed by the algorithm 
for this probluiri is shown in Sig. 3. 
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Pig. 3 Branch and Bound tree generated for the 
example proBlem 


4.3 iilgorithin - 2 

Thio algorithm is also h"''sed on branch and bound strategy 
and is siiailar to the algorithm suggested earlier. The 
branching and searching strategies used are the saiae as those 
used in the previous algorithm. The basic difference lies in 
the procedure used for finding the upper bound for Sxny node 
in the branch and bound tree, Instead of starting vath an 
initial flow of zero over all the arcs, the maximal flow in 
the problem corresponding to the nod.e from which the node 
under consideration originated is made use of. Here also 
each node of the branch and bound tree represents a unique 
assignment of arcs to different commodities. That means, 
that at each node of the branch and bound tree, the arcs of 
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the network can he c'ivided into mutually exclusive and 

collectively exhaustive sets of arcs , ... v?here 

A is the set of arcs which are not assipned to any commodity 
0 

and A^ is the set of arcs which are assigned to the ith 
commodity. Tho upT)er hound vrill he the maximal flow from 
the source nodes to the common terminal node considering 
the commodity assignments done on arcs. I’or ootaining the 
maximal flow/, wc employ the same modified labelling algo- 
rithm which vje used in the previous algorithm, the validity 
of which is proved already. 

4.3.1 Upper Bounding Strategy 

A super source is added to the netw/ork and this super 
source is connoctod to the source nodes by arcs of infinite 
capacity. For finding the upuer bound for any branch and 
bound node, the maximal flow; found out for calculating the 
upper bomd for the node from which the node under considera- 
tion is generated, is used, from this flow, the flow paths 
containing the arc selected for branching are detected. 

These paths are removed from the flow;. After this is done, 
the flow of all the commodities over arc selected for 
branching v/ill be zero. 

All the arcs connecting the super source w;ith the 

sources other than the source for the commodity to which 

the selected arc is assigned at the node ^^ider consider^ ion 

T rilff 

are removed from the network. Using the 
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nlgoritbm. In v./hich th& flow of tho ith comniodity is allov;ed 
only over the arcs in tho set A A. , v/e find the maxirial 
flov; from the super source to the comnon terminal. IText, 
the arc connecting the super source with the source for tho 
coemodity to which the candidate arc is assigned, and the 
candidate arc are removed from the network. All the other 
arcs connecting the surer source with the source nod.es '^re 
nut back in the netvjork. Keeping present flow as the starting 
flov' the maxitaal .flow from the super source to the common 
terminal on this network is found out using the modified 
maximum flow labelling algorithm. This gives the iipper 
bound for tho subproblem (node of the branch and bound tree, 
considered); 

Proof ; 

Obviously, the flow got after removing ^11 the flow 
paths which contain the canclid<?te arc is a feasible and sub- 
oiptii.ial solution to the subproblera under consideration. Since, 
in this subprobleiii, tho candidate arc is assigned to only one 
Commodity, when we apply the maximal flov' algorithm, in the 
nctv^ork in which all source nodes other than the source node 
of the commodity to which the candidate arc is assigned, are 
disconnected from tho super source, v/e got another feasible 
solution to the subproblem under consideration which contains 
the maximal flow of that commodity to which the candidate arc 
is assigned. The possible incremental flows of all other 
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comriiodities are obtained by again applying the majicinal flow 
algorithiTi on the netA'/ork containing all arcs of the original 
netA'.'ork other than the candidate arc and the arc connecting 
the super source with the source of the comnodity to which 
the candidate arc is assigned.,. Since no increnental flow 
is possible for 'iny couiAiodity , this will give the naxinal 
total flow of all coinmodities to the comrion terminal, for any 
particular assignment of arcs to different commodities . 

4,3.2 Description of the Algorithm 

The notations used here are the same as the notations 
used in the description of Algorithm - 1. 

Step - 0 j 

Add a super source S and connect each source to 
this super source with arcs of infinite capacity. Set 

1 = 0 
I = {0} 

A^ = {Till arcs in the network } 

A^ = {0} for i = i, ... p, and 
X . ^ = 0 for all ( i , 3 ) ^ ~ ...j p. 

X (j 

Solve the maximal flow problem from super source S to 
common terminal T, If the flow values arc feasible, i.c., no 
arc has more than ono commodity flovdng over it; STOP (This 
is the optimal flow). Otherwise let p^ represent the sub- 
problem. Add Pq ”^0 the List L. Go to Step - 1, 
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Step ~ 1 ; 

Check Lint L. If L = {0} , go to Step - 6. OtheTOise , 
take out the prohlaii vjhich has the largest upper bound of all 
the problems vhich are added latest to the branch and bound 
tree, let it be p^^. Update list 1, and go to Step - 2» 

St ep - 2 ; 

If U(p ) < 1 , then prune that node of the branch and 

bound tree and go to Step - 1. 0 the 2 n»'ise scan all (i,j) el, 

determine the. arc v;hich flovsmore than one commodity over it 

P 15- 

and has the niaxiinuni i x. - among all such arcs. Go to 
Stop - 3. 

Step - 3 : 

Generate new subproblcms (new nodes for branch and 
bound tree) p^^ , p,_^2’ ••* Pjnp* Subproblem, p^^^- has the arc 
seloctod for branching assigned to coramodity ’i*. Go to 
Step - 4. 

Step ~ 4 ; 

i) Solve the subproblem p^^-] ? p,^2» ••• I’rap sub- 

problems corresponding to the newly generated nodes of the 
branch and bound tree). 

ii) Trace and remove all the flow paths which contain the 

# 

arc selected for branching in the subproblem p^^^ (pushout 
routine) . 
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iii) Rtiinove the arcs connecting the svLycx source v/ith 
the source nodes nf all cciraodities other than conir.odit 3 ?- ’i'. 
Using the present flo\'.' as the initial flov;, find the rnaxiual 
flow frora the super source to ooa^aon terminal, nanmodity 

is allowed to flow only over the set of ^^rcs 

iv) Uowf put all the arcs roaio'med, tack into the network. 

Remove the arc coanecting the super source with th^ source 
node of the corauiodity *i', and the arc selected for branching 
in the subproblems Using the present flow as the initial 

flow, find the maximal flow from the super source to the 
coimiion torminal node;, considering the assignment of arcs tn 
different commodities, The maximal flow thus obtained for 
the subproblom p^^^ is used as 

Check whuthor the problem has a feasible flow. If so 
go to Step - 5. 

v) Add all the subprobloms vjith infeasible solution 

to the list L in the incre.-ising order of their corresponding 
upper bound Co to Step - 1. 

Step - 5 ; 

If f 50 t I = Go to Step - 4. 

Step “ 6 ! 

Current incumbont is the maximal flow value. Stop. 

A general flow chart for both the proposed branch and 
bound algorithm is given in Pig. 4. 
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Initialization 

I = 0 ; L = { 0 } ; 

{ .111 arcs}; == {0} 

l.= 0 for all i 

1 


k’5.= 0 for all(i,3)e ‘^■nd 
I all k 






ildd the subproblcm to List L in 

tho incroovsing order of U(p .) 

mi ' 


Ni/ 



A gonoral logic flow dingrnra for the proposed 
algorithms 


Pig. 4 






ch:j>o}sr 5 


COl'IPUTall'iOjUL ir'ERI'OK'LlZTCii; EVJJiU.lTlOJ OE 
^iLG-ORITMS :Rn]D CORGLUSIORS 

5.1 Oomputatioanl rcrforra^nce Evaluation 

The algorithms developed were coded in POETRnl'T-lO and 
imp 1 onion. t ed on DEO-1090 trae-sharing computer system. The 
programs v/crc tested over a large number of networks generated 
randomly. 

The generation of the notv;ork is done as follows ; 

Thu node-node incidence matri:€ is generated randomly. 

In the generation of the incidence matrix, it is taken care 
that the graiih corresponding to this incidence matrix is 
connected and there would bo at most one directed arc bet- 
ween any jiair of nodes. It was also taken care that all the 
source nodes nave aero indegrec and the terminal node has 
zero out-degreo. The parameters of the randomly gener-^ted 
network are the number of nodes, number of arcs and the 
nu'ibor of commodities. 

The computation'll performance of the proposed methods 
arc evaluated on the basis of the number of br^inch and 
bound nodes created and the GPU (Central Processing Unit) 
time. The maximum number of branch and bound nodes that 

2 TH • 

could be created is1+p+p +... +p'. But in most of 
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the TDranch and "bound algorithins only o froctio'i of this is 
actuolly croatod "before termination. Lesser the nunlor of 
branch and bound nodes created the efficient the procedure 
is . 

In I)1;C-1090 systen vhich is a 'lultiprograeriing 
system, the 'oroyraM ronains. in a state of execution for a 
rulativ long time i'/hilc the central processing unit 
takes it UT 1 only in > ermittnatly . 

I'or each co-ibination of the parameters of the network 
five arobloMS were generated end programs for algorithm 1 
•uid ••‘igorithn 2 vfer.; tested on theso problems. Problems upto 
the siKo of 200 nodes 300 arcs and 5 commoditios v;ore solved. 
The comput<ational experience is presented in Tables 1 to 4. 

Th-.' nu.'iber of fc-isible solutions to tho problem token 
uc iii this, ■^'/ork is proportional to p whore p is the number 
of CO iiodities 'Pil ■'.i is the nu ;ber of arcs. So ns the 
number of arcs and the nuiiber of commodities increase, the 
computational complexity also increase, lypica.lly in n. 
branch and bound procedure only a fraction of the feasible 
snluticns are crmrierated. In the suggested algorithms, 
because of the tight upper bound, the number of feasible 
solutions enumerated explicitly is expected to be very 
small and so moderately sized problems could be worked out 
using these algorithms. Most of tho problems tested were 
solved in less than 5 seconds of O.P.U. time. 
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Algorithm - 1 v;as found to be computet ionallj suporior 
to Algorithm - 2 as it has lose memory requirements and also 
t'-lecF leas C.P.L). time. In almost all problems Algorithm - 2 
tooic much more time than Algorithm - 1, But it was noticed 
that as the density of the graph increased, the rlifforence 
in C.P.IT. time botwocn Algorithm - 1 and AlgorithPi - 2 
d G creased . 

5.2 Oon elusions 

Per the problem of finding the maximal arc disjoint 
flow in a multicommodity network considered in this vjork, 
two algorithms are proposed. Both the algorithms are based 
on th e branch and bound philosophy. the upper bounding 

strategy employed in both the algorithms violates only one 
constraint (i.c. the arc disjoint constraint) of the ori- 
ginal problem, tight upper bound is obtained. Phis explains 
for the low nuiiibv-r of nodes created ixU the branch and bound 
trees generated by these algorithms and their high computa- 
tional efficiency. These conclusions are based on the 
experimentation conducted on a largo number of randomly 
generated networks. In aliiost all the networks generated 
for experimentation, Algorithm - 1 was found to be computa- 
t-ionally pior^’ cfficiait as it took less C.P.U. time and 
required less core moniory than Algorithm - 2, 



45 


5.5 Scope for future Work 

The oolutioa jroceduros propoced for the clans of 
disjoint flows considered in this v/ork are hased on the 
hranch and hound technique and so implicit cnuineration of 
all the feasible solutions -'•re done. It may he worthv/hile 
to think in torris of exploiting the comhinatorial structure 
of the problem and come up vdth more elagent methods like 
labelling method. 

The present work was concerned with pure disjoint 
flows, i.e., vdicro all arcs or nodes allov; only ono commo- 
dity to flov? over it. further work could be done by con- 
sidering the mixed disjoint floivs, i.e., where some arcs or 
nodes allov; simultaneous flow of commodities. 
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SHRCE 


STARK J) 
FTNTSP(‘T) 
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f^'ut’Per of Modes 

«u'^ber of Arcs 

:'ii'u’’ir)er of Com«iOdities 

Nodenun’oer of toe super source 
(tols Is eauaX to •''fJOc.F) 

Rtartinq node of tfe J th arc 
Flnisiilna node of toe J fh arc 
Capacity of the j fh arc 
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- The source node ot the 1 th ronifodttv is 

numoered ' t ' 

- The arc connecting the soper source and the 

source of the 1 th copoodify is nu'^nereu 
*J' and tne capacity of this arc is set 
to be a larae numberfinflnlty) 
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T Oils 2 ? 

riPFh ‘aiNTTsIlN, Frr,K='NRTGEN.inu'3^ 

OFFN (UNlTsinU, PlhE s'NEtnt 1. .DAT') 

REAUfllN,4)NMU0K, 4ARC, NCOilM ,SINK 
READdTNKlCSTARTC.T) jJsl.WARD^ 

READ C 3 T » f H tl nf '1 ^ ’ 

READCXINK) = 

DO 3 I»l,NARn 
TFLOWCDsO^ , , 

DO 3 Jsl.NCOMM 
FLOW(J,n = 0 

TOTCANaO 

SORCKsNNnOE 

This part Of the program .using the Inout data 
generates the following variables which provide 
a better data struture to the program 
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’■$ 


NeTCK) - K th element of the vector of 

nodes NET 

degree (I ) 


The total degree -of each. node | 
<le, total no, of Incowming R 

arcs+outgoing arcs> 


CLOECi) 
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• The starting point. in the vctro 
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s assist 
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0 
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0 
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PDC&. fK) 


■ net of nodes whicn are dtroctlva 

• connected to tne node 7 o 

■ The edoe correspond.incj to the 
K th element of 

■The direction of tne edae 
corresnonolno to the fs th 


element of hFT 
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C 

c 

c 


c 

c 

bO 

SI 


c 

53 

C 

52 


F, = 1 

CLUEt n = 1 

no 30 r S 1 ,'‘dnDE 

m 20 vT = J,NARC 

IF f start GO Til 10 
IFfFlNTSR(J3 .K0,T )Gu on S 
GO TO 20 

5 dET(K )=STAhT(d) 

niRECT(K'»=2 

FDGE(K)sJ 

Ks*f + J 

CO fO 20 

10 )=Ft <TSH(J) 

nmecTf Ki = i 

FDGEf Kls.T 
K=K + 1 

20 Cuoi'tw'iE 

CL(lE(l + n=!< 

30 nEGKEE(I)sCLOE{ i + n-CLOEf 1 ) 

TNCUMBsO 
INDCAN=1 
DO 40 .J=1,NAPC 
ASSIGN (J>=0 

IF THE MAXIMOM POSSIBLE FLOW IS FEASThT.fi (ARC-D 1 SdOlM-n 
ALGORITHM TERMINATES AS THE Fl.CM TS THE MAXI'^AL ARf- 
DISJOINT FLOW 
JAsO 

CALL bound (UB.K) 

IF {K,E0.0)Gn T0998 

IF INFEASIBLF ADDED TO THE LIST, 'K* IS THE ARC SELECTED 
FOR BHANCHTNG. 
nn m .tsi 

LSTCANdNDCAM, DaASSIGNCJ) 

LSTCANClMDCAN,rlARCtl)=K 
LSTC AN ( I N ()C A M , d A RC + 2 ) = tl H 
TOTCAoisTOTCA^' + l 

IF THE LIST TS EMPTi ALGORITHM TERMINATES 
IFf lNDCAH,Et?.0) GO TO R9R 

IF LIST IS NOT EMPT’/ THE l.AgT PROBLEM In THE LIST TS, TAKEN 

DO 52 Jal.NARC 

assign CJ)sLSTCAN(INL)C AN, J) 

IF THE UPPER BOUND IS LESS THAN THE INCUMBENT 

that node is pruned. 

IF(LSTCANCTNOCAN,NARCt2).LT.INCUM0T GO TO 54 
11=0 

CANAOITE PPOPLEMS APE GENERATED HI ALLOWING EACH OF THE. 
DIFFFRENT COMMODITIES TO FLOW OVER THE ARC 'K' 
KIaLSTCAN(IN0CAN,NARC+1 ) 

INDCAN«1N0CAN-1 
DO 55 IPsl,NCOMM 
ASSIGN(KT)?*IP 

DO 49 NsI.NCOMM ^ ^ 

^ t 2 


55 ^ 



no 49 Ms!l,MA«C 
49 PLOWCN,M)sO 

CALL BOU^»a(UB,K) 

IF (K.Ms,o) no TO m 

C IF THfc FLOW IvS fPASlBLK Ti4CUMHFivT fJPOATFO i *■" THF 
C UPPER BOUWn TS greater frU;*' the TKriO*Br'.T,ni-FP;P>/TvSF 
IF (UB.GT. T.mCUMB) TMCUMBsUB 
IF (f IP,P’Q..<CO»^M).'\hD.(IT.Fg,0)) GO TO S'? 

GO TO 55 

C TsvFiiASlBLF. PRDPUEi-iR ARE AijOEP TiJ THF 1,13^ T-^ T'L*. 

C nPDFP OF THk’ ijppFO BOU’Cij 

57 TI=II+1 

no 59 Js1,MaRC 

59 I ASS T G C I T , 0 ) s aSKT C 0 ) 

SELAPCf 1I)SK 
UBOUNDdT )=tiR 

55 CONTINUE 

no 56 1=1,11 

56 IR(n = T 

C THE INS!, Poni T.^P VSRT*^ '- IS UtoKO V' fl ,vG T:, '.i Uh_ HRliHOPT-wS 

C the INCHFABi ORDER uF I’itF UpPE» HuU"H> 

CAT.!, VSH'^P * (!jPOH*D,Tf 

no 5R Jsi,!! 
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IF (TNnCAf-..GT,50)G!l TO 617 
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DO 60 Ia1,NARC 

60 LSTCAN( iMDCAM, I)=IASSIG(THCJ) ,T) 

I.STC A N ( 1 H DC A M , N AKC + 1) =SF.L AFC ( I P ( 0 3 ) 

58 LSTCANClNOCAM,MARC + 23 = UBnuMD(J'> 

GO TO 5i 

54 TNDCANaINDCAM-1 
GO TO 53 
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999 WODENsMNODE-I 
lARCNsNARC-NCOMrt 
CALL RTXWE(TIME2) 

MSEC=TIME2-TTMEl 

ARCNslARCN 

ODENsNODFN 

DENSTYaARCN/ (QnENl= (ODEW-1 3 3 

WRTTR(TQU,794 3MiinFTI .IvOCEN, I ARCN,HCOVM,IUCUNR,TfjTCAN,f'SEC,DENST¥ 
794 FnP.4ATC/5X,7T8,5X,F6.3) 

618 CLPSE CUNITsTI.')) 

CLOSE CUNITsTlIU) 

STOP 

END 

C 

C TOls routine finds the nsax.flow and check for the feaslollttv 
c of the flow . 

C FLnw'(K,I3 * Flow of the Kth comrriodlltv over the 1 th arc 

c U8 " Maximal flow 

c K - The arc selected for branchinq<ie, ,of the arcs 

c which allow more than one commodity to flow, 

c . the arc which is havlna the larqest flow> 

c 


CUB.K) 

?A-Z) 

NMODE, N ARC, MCOWMiSORCE, SINK, START (300), FINISH (300)* 
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SUBROUTINE BOUND 
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common ' 
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" THE VECTOH of ^'ODES WiFTCU AilK 
-NUHBFR OP’ NnOES .vHtCH APPj LAhFlJ.EO 
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1 Tr ARC IS FURS'- ARD 
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n«? tse i th ;)UijE luniCATt'-’G i'‘?p: arc ay '■.■htc! 

'lijm-. i tjni' ?,A;iEi,!,,b;n 
L5CI) -uAbFL UP’ TIE I TH i THF rs-jDICATlEv, tHp; Cn^f'.jniTV 

I.ARGF«R999 

fIB=0 

no 7 i=i,'!Akc 

TFT.OWCT )=0 

no 10 Lsl. ,NC0 1M 

T.KliTsT.ARGF 

L2CLT=SURCF 

L3(LT=1 

T.4(L)aT, 

Ii5Cl.)=r. 

lilSKOTsL 

continue 

NLTSTsNCOMM 

ALL THE SOURCE NODES ARE INCLUDED IN IRPT LIST OF 

LABELLED NODES AND THE SUPER SOURCE TS HOTU 

LABELLED AND SCANNED 

DO 14 laMCOMMf 1 ,NNnDE-l 

L1CI)*0 

L2Ci)«0 

L3CITS0 

L4(I)=0 

L5 ( I ) aO 

USTCDaO 

CONTINUE 

COUNTsO 

THE NODE NriTCU NAS FIRST LABp]I.LKU BUT STlT.Ii UNSCANNFU IS 
FOR FURTHER LABELLING 
KKaLTST(C0".'.T + 1 ) 

TFCKK.FQ.O) CU TO tOO 

IF NO FURTHER T,ARELLlsJG IS POSSIBLE MAX FLOW IS REACUEn 
JJaCl.UF(EK) 

LLaCLUECKKl+DBOBEeCFK)-! 
no 50 MMaJJiM, 

IlaNETCMM) 

IF CHCII).GT.O) GO TO 50 
NNaEOGECMM) 

rF(CA5SlGN(NNl.NE.0).AND,(ASSIGN(UN),NF.L5(KK))) GO TU 50 
IF CCDIRECT(MM5.EO.l).AND,CCAPC!4N).GT.TFIaJH(NN)))GO *’'0 30 
IF ((DIRECT(MH).EQ.2).AND.( rFL0MCNNj.GT.0)3G0 TO 40 
GO TO 50 

' 14 ' 


CHousp; 



C IF FURTHER b'^UFLLING IS MOT PnSSTtT-F J»'ST nhCftHs^' 

C THE ASSIGnMEMT OP THE ARC IS I) fFFt.REMT, ARC 
C STORED SO THAT RATER THE COMMODITY TO WRTC'? ThAT APC 

c IS assigned to cnubo he trifh with 

C T/ARELRING OSING FORWARD ARCS IF aDDiTrnr.tAL CARACJiTE IS 
C AYARIHLE 

30 D I F pssC A P ( )'J H ) - J’F LO W ( N H ) 
f.l(iT3aM IKK) 

TFfDTFF.M.bl (KK) J bUIDsDlFF 

Ti2r I.T)=KK 

H3(IT)sl 

T,.4(iT)=WH 

r,5CII)=L5tKK) 

Go TO 45 

C t.ARKIibTG USING oACKi-ARn ARCS IF VhO.> OVRP THE ARC TS >D 
40 I.) fIT)sLl (PK) 

IF CTFr.Dw{H!O.Oi’.t.,UKr,)) LI 1 1 1) sTrirU: ( f-Li ) 

L2CII)sKK 

L3(IT)s2 

I.4(II3sNN 

L5(II)=L5CKK) 

45 IFCLI (SINK) .GT.COOn TO 55 

C IF SINK TS LABEl.t.EO FLO.v IS AUGlH'LNTFn RY AriOU^'T Lt (SINK) 
no 47 T=1,YLIST 
TFtJLT.FO.LTSrCT)) GO TO 50 
47 CUMTlNWe 

NLTSTsNLISTtl 
MSTCHMST) = TI 

50 cdntinwe 

COUWTaCOflHT + l 
GO TO IS 

55 CHANGE*L1 (SINK) 

UBsUB+CHAnGE 

L=S1NK 

60 IFCL,EO.SORCE)GO to 12 

MaL 4 a) ' 

NsLSCL) 

IFCL3Ct.).E0.2)GO T065 
C FOR FORWARD ARCS THB! FLOW TS 
C ARCS THE FLOW IS DECREASED 
FOOWCN,M)=FLOw(N,M)+CHANGE 
TFLOW(M)=TFLnw(H)+CMANGE 
IjSL2{1;) 

GO TO 60 

65 PLOW f H f M ) sFLOh ( N , i^) -change 

TFT,OW(M)=TFLnwf.M ) -change 
0*1.2 (Ij) 

GO TO 60 

100 IF ( JA,t:o,Ncni4FO on to 70 

JA=JA+t 
NOISTsl 
I.ISTCDsJA 
DO 11 TJs2,nC 0'1M 
U T.ISTCIJ)sO 
rn TO 1 3 

THIS PORTION OF THE PROGRAM IS FOR SFLFCTIMG THE ARC FOR 


INCREASED AND FOR BACKrtARD 


C BRANCHING 
70 


MEDIATbO 
K *0 

no 80 MsNCOMM+l ,NARC 
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TFfTrL0wf’O,p:a.0) th so 
DO 90 N = 1,MCDviM 

IF (rLni^fN,Mi,FO.TFr/n»(M)) go to ho 

TFfMFDIAT.Gr.TFbn,-J(H))Gft TO 80 
90 CDMTTNtlE 

MEniAT=TPf.nftf*vil 
Ks M 

80 cfjNTrruHh' 

RFTljWN 

Fi\iD 


1 6 



UUOUUOU’ 




S ft: ^ Si 


Ai..(ir,RX fsj-. 


& Si /v S' ^ R 


•^►•.''' 'v: >•!..'!!!;■ rRi'jCf-^n’-^Pi,.; mh; .sni.^'viv; A 

' ~iv:r!4f. ARC-':'i.sjni.'i r ■uj.i, r^^l..■^■■tf^/l j.y 

i'!«: ^I. I 'lAi;, KuO ; PPOf’tf*'' 


T h f.i.iA'O )A J 

n n(n } 


T 

i'Tl.lCT 1 

I 

'i i 

a" 

GE 

fk 

(A- 

Z) 

C 

n,-' V 

, j , 


* 

•V A 


' 5' 

•li-/l,.SuPC 

1 


;^(V 

":i 

JP 

OF 

tjP 

d] m 

f r.'h:,r.r?o 

T 

f. 



T» 

(3 


) ,1' 

1; J.hu(3i.'; 

1 

i/KCUEFC 


0 ) 

f 

Cl. 

(J 

d') 

0 1 , FDGF ( 

7: 

50i.CU, 

f* i f 

‘ '' ’‘i 

e 

31 

. 1 K 

f n- 

u 1 .1 ■ ■ i 3 U 1 1 ) 

3 

,(.STCA.'i 




,3 

\i 

) / 3 

"h.'.i-vOC? ) 

4 

, TFLOFt 



n 

) f 

A,S 


(, 3u,i ) 

RfeiAh DFi4 

STV, 

n 

0F> 

.• 1 

.•UaC 

ii 


ft 4 



II If 


^ ff 

!< ,1 1 B b; (1 1 (I 

# 

f'-h'uDE 






Mi 

fU:er of. 


i.'AhC 

uCu'^fM 

sorce; 


# STARIfJ) 

« FTNISHCJ) 

# CAP(J) 

# ^fOTfe; 


" Giriner of Arcs 

- ^!u‘Tic>er of Cor, uvo(j^ ties 

- ^'odenumoer of tr.e surer source 

(this Is equal to uuonp) 

- Starting node ot the i to arc 

- FinJsnlnq node of tne vl tn arc 
" Capacity of the J t-n arc 


i - The source node ot the i th romonltv Is 

I numbered 'I' » 

# - The arc connectinq the stiper source and the n 

» source of the i tn cott.mod.ity is numperedS 

# 'J' aol the capacity of t,hts arc is set # 

# to be 3 large number (inf Inlty) « 

Tlds?i 

TUna^2 

nhFh (U uT ,'SM 1, FU..G=' /U'l’f’hF)- .mu' ) 
opFu lUMismu, yjhii =' ■i!''tbib.r)AT') 

TuTunhsu 

REA V ( II it. , f ') - < h E , ■ ' A isC , Y C 0 K >. , .'J I ; J N 

PEAnfXTW, + KRTABTC,,T),Jsi ,,JAHC) 

RBADdTN,* ) CFlII^lrlCJ) ,.l=l ,NAnC) 

RtADCIIN,*) tCAP(.J) ,jsl ,hAKCJ 

DO 3 lalfhSHC 

TFL0'4(TJ=U 

DO 3 Jal,hCO«M 

FLOWCJ.IIsO 

LARGE5S9999 

CALL RTIMECXIMFI) 


This part of the prograin msinq the Inout data 0 
generates the following variables which orovide 0 
a better data struture to the program | 

ffel'C |iif )‘ "" ""f t"h" "element" ' 

nodes NET 

DEGREE(l) - The total degree of each node 0 

'"2 i ' ^ 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


la 

la 

la 

ra 

a 

(a 


f \ 1 

k^.) 


<ie»iotfii r^o* 

a r c s 4 o 1 1 1 a 0 i v ! a r r s > "a 

-* start in<i potT^t tn f?ip vctr-^ 
ot p.oti^^s whici'i are olreetiu*^ 
CQpnocteH to t^**e no^ie f i 

'■»*oe Oiine c<^rresona^Lt nq to t:r»o a 
K tn eie^ie’^t of 

-•Tne qJrectiOii ot tne edoe *a 

corresi^onuioci to iho, to a 


ejo;>'.eitt' ot Ff 




Ks:l 
CLf^Ef, l)s1 

nu 3f> 1 ::‘'t,''- 
raj 2 ^ J r 1 ■ 


/n)K 




If. ^>^TA rt( a 1 i ) ■’"u 

TF fFTtyTSHCJJ.yn.r hio ’! 
r,0 TO 20 
5,„ « Nti'^'CK:) = &THfn(il) 

0iRFCl'(K)s2 
FDGE(K)s:J 
KsK+l 
r,0 TO 20 

KDj,ErK,=,J 

20 , CONriNi'E 

CliUE( I + l )SK 

^ hEGBEe(i)aCLIfE;(l-H)-CLUF:f J) 

INCUMsO 

INOCANsl 

UBay 

MAXIFOM POSSIBLE FLOW IS FEASIBLE CRHC-0.lSj,niOT) 
_ _ ... 


'Kl' IS THE r^nC SFuFCTEO 


C algorithm Tfe'RhJNfli’ES AS THE FLOW 
C niSJOIfIT FLOW 
TSFTay 

no 34 Jal,f!ARC 

34 ASSXGHfJ)=0 
OAaU 

C A L L , B 0 { J H u ( U B , T » I S F: T , A S S X G 1 0 
IF CTK.EO.f'I GO r.j 99 
TOTiiODsiT) TMOni + l 

S I**’ I”t'’FA.SiBLS ALOED TU T-fE LIST. 

C FOR, r^RAuC {TAiG. 
no iS Ja1,f?ARC 
LSTCANC £W[)CAM,6» JlsASSKOHJ) 
no 35 Kal-MCOMA 

35 LSTCAH(IMOCAM,K,.’f)aFLO>'(K, JJ 
LlSTClNDGANfl 5=1K 
f.ISTClNDCAN,2)sUB 

LIST TS EMPTY ALGOHITOM TFHMIMaTKS 
40 IFCiNDCAN.EO.OiGn TQ 999 

C IF LIST IS MOT EMPTY THe: LAST PROBLEM lA THE: LIST IS TAKEN 
DO 45 Jat,NARC 
ASSIGN(J)8 LSTCAnC1MDCAN,6,J) 

TFLOWCJJaO 

DO 45 Kal ,NCaHM 

FLOWCK,J)sLSTCANCIMDCAN,K,J) 

.'22 



c 

c 


c 

c 


45 TFtjOWCvl)=Tf''’L->v,f J) 

IF (r,isifi‘4DCA'^:;j.f,F. r.jcur ) m to s"; 

KIsLlSTdonCA.'^a) 

(JB=LTSf { r one , 2 ) - TFr.n 0 c x i ) 
i:!gDCao'lsI''!()rAM-1 

ALL THF 'J'Lnft PATH?? .HlCi' CL. TA T.:-4 T^IF nf’C fJr"'. t 'L-'"' 

FRn,-.i FF.nw BY CPF !."liTVi;K 

CALI, pnSM.HKKl,) 

Tonf,:x=0 

no fc;B T=:1,'.‘*A»C 
TTFT.n^f lIsTFr.frn T ) 

DO foS IKFsl 

65 1FLf:r”( TKT , T J =FLO ' C r k i , r J 

Iv'DsDft 

CA^'JAOJTF pDDBi,,>v,$ 'Lv ■•SR ,v'T;n BV .'J.iOL*. T ■ r, i sV ?■,*■ 

DlFp'FHFrAr CU'^'-iOUirTF^J 'r/i 

DO 55 Tf> = l ^^vCCr-".! . - . 

AagT'L' f KT J5S,IP 
flBsilTB 

DU 77 I =1,0 ARC 
TFL(l>»(T)siTFLO%l I 5 
Ou 7? JFT = ,1 

72 FLio, ( 1 K i , r ) =■» Ff-tjo t Fr, r , r ) 

CAPKIsCAPtKI) 

I6FT=1 
Katf 
•I A"0 

CALL BOU*IDCU0,K,ISFT,A.S6Ii^B) 

IK = IP 
CAPlKDsO 
ISFTS2 
JAsft 

CALL BO(INOCUa,lK, I SFT, ASSIGw ) 

CAPCKDsCAPKT 

IF THE FLOW IS FEASIBLE IDCUMBflnT IS Ut’OAiEU THE 
UPPER BOUiilD IS GREATER IHAM THE TtNr’J*O.F,.iT,ni'1.|:.RdSF PH'L'^EO 
IF CTK.MF.nj GO ru 75 
IF (UB.GT.TNCUM) tVCUn=uH 

TF((TP,EO.rlCOHM).A^n).( roDEX.Ey.o)) OH ,*0 
GO TO 5S 

C IREEASTBLE PRUPLFmS ARE ADDED To THE I IST To THE P.'CPEAST liG 
C OKDER of the upper BOlIfNlU 
75 TNDE)f=TNDE3F + 1 
DO BO -.Tisl^MAHC 
T ASS TG f I HDEX , 6 , U ) s ASS 1 (L* ( J ) 

DO BD K si, neon;-' 

BO TASSTGC lMUFXfK,U)=FLO/-(fV,in 
SELARCf IDf>FX)sTX 

imOiiDDf lf?OEX)sUrt 

55 CUTTIfJlTE 

DU 85 Ts1,I«DEX • 

H5 TH(1)=I 

C THE I'HSL RUHTTf'iE 7SPTPD IS USED ARRALUTnG Ih T'HE PRUfluF-iS T 
C THE TWCREASIMG ORDER OF PHE UPPER RDHOD 
CALL V5RTPM(rJBnU0DrIWUEX,IR) 

DO 90 JBI, INDEX 
INDCANsIMDCAM+l 
IFCINDCAN.GT.SO) GO TO 617 

totnod=totno5+i 

DO 95 fsl^TARC 


C 

C 



u ouuououuuuuooou 


liSTCAMCIMDr'AM,6,T)sIASSlGCIlUU) 
DU 95 K=1,WCnHW 

95 I.STC N C I M DC A V , K , U = I 435 1 G ( T k ( 0 ) 
T. I ST f X woe A W » 1. 3 ® vSK u A G ( 1 1- ( a } ) 

90 i»iST(:iwoeA’‘^?j=a« jukDCj) 

GO 'xn 4 0 

50 Twr'CAf'isiMnCAV-l 
GO TO 40 
99H T^ru 

999 MoOSOs'k'/Me-l 

J APC"’ = Mh'JC-'‘ " fi-'''- 
CA.i.b HT1 "h;( n.’.>!21 

ftKC-'-sIArtC'.: 


# 5 , 1 ) 
# # T ) 


794 
&1 H 


OKwsT'vsAPc^' / Cijr*;'’'' f i ) ) 

o R T -J, h: ( T 0 r > , 7 9 4 T'n ■ < ' p T ; .. O D K ‘ , 1 A R G . i , , 

A T C /R < , 7 T , 9 V , Fb , ) 1 
rLn.SFilIfkC J'sj .1 O 


ro'b' , 


"i" I iF' O' 


■)F 


ci.n5P:cnNTT=iou) 

STOP 

Ewn 

This routine fJnHs the ii .*» x , 1 1 o =v anr! cnecF for the f ea;,1 f x 1 1'' 


,'Y 


C 

C 

c 

c 

c 


of the flow 

FOnwfK,!) ~ Plov. of tue Ktt coFFOdltv over the r th arc 

HB - Maxima) flow 

K - The arc selected tor L'ranchina<ie, ,af toe arcs 

Which allow more than one roouodltv to flow, 
the arc which is havino toe larijost flow> 


SUBHOUTIWE BhUNI) ( (rB,K, ISbT, ASKlGi ) 

IMPLICIT IWITGFR (6-Z) 

COMMON w W on K , W A kC , C n >; M , 3 URC K , 3 1 h t\ , 31 A P X , F T (>' 1 .;Vf , ! 
i ,WET,DlPgCf #OFGREF,CLn£,onGF, TFl.n ,M.’' 


'Ah 


l /■'IfcT.U 

INTEGER 



'' hC , 5 I'A r-;T ( 1 .i n ) , F f 1 3H i 1 ''(U , 
Cfi'iF C7>)h) , 


LIST 
WLTST 
COhkT 
Tdf 1) 


L2ri) 


T 3 (l) 


L4(I) 

hmi 


- T),e vector UF i"i,nES kdlO? AHF l-APEt-LFp . 

- lOMBFR OF i.ODFS VHlCti ARf. l.AbFD.hh 
-wtlvBFM OF ipuES -OnCH flPF A'-'h SCA.jtED 

-l.AdFi. i)F ThK T Th -ODE 1 FPirATihG XHE 
fF MOGf-iE.d'ATXri;. pOdSrPbF 
-l.Ai'Fi. triH TPE i Th rbjhfc. r,Fl>Tt7 rJ^'G the MJOE FRUM 
/HiCk T'i GOT f.,.Ai'.FLTKn 

-i,AGEL FOH toe X TH YUhfc. IMpTC AfTr'iG THE niRECTTOW 
OF THK AJ'C kV which KfjhE i GOT L.AkFl.t..En 

1 IF ARC I>S FOPFARP 

2 IF ARC IS RACHWAPL' 

-LAHFL of the T th kDDF IbhlCATlNG THE mRC HY '’'HICR 
WORK 1 GOT l.ABFLLEh 

"liBkSk 0*^ THE I Th RODE THhlCATlOC, T'-'E COhMutHTY 


LARGES9999 
12 DO 10 l.=.1,WC0HM 
l.l(L)s).ARGF 
L2(L)=S0RCF 
L 3 C L ) * 1 


7 4 



T,4fL)sTi 

r 4 S(L)=t 
10 CONTTNUe 

NLTSTssMCn^M 

»i 4 f. THP SOnBHfc 


C 

c 

13 


.nf)r 


!».» >4 


PE 


I.ABETjbKa T|Ui; 

T.Awrr,{..Elj) Ai;u rtCi. 

Du H TsMC''‘' -* + t 
t.lf ,t)=0 


i.f T. 

b^iiFCK 


■i'”E '■’(> 

I,' 


inoK' 


1 2 ( Usf) 

I. 3 f ). IsO 
h 4 (l ) =0 
T.fj f i ) =0 
T.ISTf n=ii 
14 CU'^iTr'^nE 
CUfft-.Tsn 
IFCl.S 6 :T.Ef 


I 't 


'‘fV 


IFf ISET.f'w.^li.; ) 




i^AX j H Pi'.Av.J-'p,''' 


1 1 ; 9 
I '**' t 

ni; n Ls1,rClMM 
11 i.isiaisf. 

C THE Nnor* 'JHicn ■'■hS Tu^’':ET..ii-i' «"r 

C FfJKfHEk 

15 KK = LTiiTCr:uH .T+U 

If CF^.EO.oi «n ra jof? 

C Tf HH FUPTMER uG .1 « VTuS.Sirn.F 

1« JJsCI.UF(Kk) 

i.i4»CT.uE(KE)+r>E«HEecisK:)~i 
DO 50 Hi4sjJ,0b 
Ils.'iEKMMJ 

IF (Mfin.GT.O) :;n tu so 
NNsEDGECMm) 

ifccap(ai''?).eo.o) go ru so 

IF tAS5IGNCN.^’J,ME,O.A)MD,ASS10MCriO),?<.f--.T,5fK'':)'> 00 TO 5u 
IF (fD.lRFCT(viM),EQ,n,Ai'-D.CCAHCi.'''J.GT.'’>T.<'i-fl ‘ .<0 

IF {(DIRECTtMM).E\^,2),Al-.n,CTF.bn,.!(’ Tj 4P 

GO TO 50 

C iABEOLlNG HSTNG FaRWAR!) ARC;" IF AODl T 1 0,: A*, CAkACITTt IS 
C AVALIBDE 

30 DIFFaCAFC:|fn«TFL04(f4 0) 

I.lCIT)=Lt IKKl 

TFf0IFF.f/f ,01 (Kkl) LI Cn)=?nFF 
I.2(IT)=Kff. 

L3C1TJS1 

L4(;n 

LSn DsLSCKhl 

G(f TO 4b 

C LABELLIG USl»’G t;AG»<,sAKO AkCS IF FlO'.: I'VED TtiE AkC T,s >0 
40 T.mDeLl (Kr. ) 

IF (TFIO'UOM -FX.Ot (Kk) ) 

i.2(it)=k:k 

L3(IT3s2 

l.4CII)ttNW 

L5ClT)aL5CFKI , , 

45 TFCKK,Ei^.snRGE)L5Cn.) = Tl . „ . 

C IF SINK IS LABFLLED FOn.' IS AUGUMENTFO QY AmOU'-.-T LUSIWiO, 
IFILKSIMKI.GT.OJGD Tt) 55 
DO ^7 1=1 ,NLTST 
IFCII.EQ.LTSTCI)) GO TO 50 
47 CONTIMTE 

WLISTaMLTST+1 ^ ^ 

2 5 



LISTfNMSl ) = IX 
50 COHl’INUK 
57 CO?m'>’=COtiN'^+} 

GO TO 15 

55 CHANGEkOI CSIM/O 
npaliRtctfAMGfc; 

L=51«K 

60 H*Cli.E0.5a??CF) >‘Li i> 

^=T,4f 0) 

’’=15 (]..■» 

C FOP FOPv.AiOi (in'CS rrfc. HT.r:;.' f:<; I'.CPfc Ai-iKi.' n' )•' fOH tif- ;; X ^ 

C ARCS THE Ft.Oi*' rs ;jKCr:EAEF(5 
It (1,^(1..) . f, ro 
FLOft ( , .> ) =Fi'.p'.., f :■ , ■ . ) 4Cr ^ ■ OF 
TFOO M t '" ) = VF L 0 V ( M ) f C >- \ G F 
T.=G2(i.I 
r,ij TO 6i> 

65 Fi.n,- f , , ■iTsFv.'O., {' ,1 

TFTiUF C ) ='l F 1 . -1 1 -G; . A G?i, 

I=I.2(L) 

GU TO 60 

101 ''.LTvSTb*jco-'.‘1-1, 

LL=0 

Dfi 103 L=i , 'fCu'G'i 
TF (I-.FQ.Kl GO 10 103 
T.i.b=LTi+ 1 
IilSTtliLI=0 
103 COMTTAIUF 

GO xn js 

102 Vi^iiiTsi 

LISTC1)=F 
GO TO 15 

C THIS PORTION UF THF PFOGHAF T6 FOR SFLFC'’MGMT IRK APC Kr’R 
C BRANCHING 

100 TFCJA.Fg.NCuMMUO TO 70 

JAaJA+1 
NLlSTal 
LISTfDsJA 
DO 111 IJ=2,NC0iM 
111 t|S.I(II) = 0 

Go TO 13 
70 MfcDlAIsO 
K=n 


DO ao A»C 

vi .Fg.u) GO rn ao 
Dfi ‘ID M = 1 , :;CD G'’ 

TF c FOOW r, Kb, I'FGD i(s)) GO “nJ H!> 

XFrMKUTAT.Gr, i'Fj., :T )Gij ff* P*) 

90 COWITNUK 

MEDlAT=lPhn.'; ( 1 ) 

KsM 

80 CONTIKHK 

■RftTuR-I 

END 

C %%%%%%%%!; 

C THE PODTI'.rF PMShOU HF;.'nVKS AGO THfe; KT.DK DATH5 WHICH 

C THF arc '4' 

SUBROUTINE FnaHUD(K) 

common NNnDE,NARC,NCnMM,SURCE,SINK,STAPX,FIWTSH,CAP 

2 6 



1 ,NET,niRKCi ,i)KG'^ci;R,Ci.''if;,eOGK, 

INTEGER S'i'AKT (3U01 IT’-i { 3uO ) .C « ;iO'' 1 , ^ r ( .'i-rO j , fiui'C'n , 
lDEGREEC?nn) ,>jr.,u,v;(J.oo J ,tnGH:(6<^A) 

2 s n RC K , N c u fr* St \ K,nL a - ( 3 o ) , f f . i . ’n s , 3 v.- r j , ( :•, : ' : j c - > ; 3 •■' ' / ) , c ' g, /, ^ 

00 5 KCO'^si f^'cn.R' 

IFCrT-.OWCKCf^'^fK). FO.cn GO xc; “S 

CJi 

I..|WsF[/lwf aC.O ,K) 

TRFOTU 1 )=•'. 

??■’? .0 w I K c n = I F r, a - ( k c ,¥■) 

F ii I'M, f K c 0 C *' i) s «' i., . c f . c o ■ , K c n 5 V , f K. r , I " , r-n 
TFCvHTARrtn. ni.KC.fn'JL- r^* io 
C H a ..t GE = C 3 ■■ ; c; + *•• « f v. ( -■, c . X- , '< ) 

GU TO 4'^ 

10 TA’V.VsLTM-C'O: '-GF 
TI = 1 

Nf;sK 

15 F,K=m'AKxf -n 

JJ=CLUF(fk) 

l,L=CT(UF( FK 1 tDEG^FECh 4 jt-j 
DO 20 vv,s!,,f,i,r. 0 
a?k=En&FCo n 

IFfDTRFCTC^M- n.RO.l U,n Tl. 2C 
DU Hi TIKs1,IT 

B1 IFCNM.Eg, iREDurilK) ) on Tu 20 , , , 

I F f t F L 0 .-4 f !< C U N ) . ! ' F , cn , A I . D , f F I I ft Ft ^ . F, * ■ . 0 T A R i' f I E ! ■ i ' H, I ) m 

1 GO TO 30 
20 rUNTINUt 

30 I F- ( F 1. 0 v< ( K C n !-i , ('J ^n . L T . I A T ) I A ^• T = f j . f • v* f K 0 g . . C 

TFCIAMT.r.E.O) GD TO 4f) 

IIslI+1 

IREDlfCIDsWiv 

IF<STARTfrt''n.wF,KCrF-l) GO TO 15 
42 CHANGE=CHARGF+rAMl 
DO 50 Ts7 IJ 

TFLOW ( IHFDU t T ) ) sTFLU w ( J BFOn I T ) ) - 1 A ^'T 

50 FLnw(KCO'^,THFuntI)} = FLOii(KCO'' nf^Ff'KT n-T Vf 

IFfCHAMGF.Di'.blM) Gl,. TO 1 '> 

40 FLOWfKCOM,r)=LTM 

CHAFGFsO 

TFfFTNTSRtHn ."jF.ftn'A) Gg TO 41 
CHA''jGE=I.T,'.. 

Cl, I ‘i'f’ ft 

Tl IART5^^JT/l-’t:0.’\^ GP 

Tlssj, 

Wv;sK 

45 FKrFrin;o>-’(' .n 

JJkCTUF C-'O 

T.LsrCLUF C K !•. ) + DOGFFt; ( f- K ) - I, 

DO 3*^ I, 

■nnseogfcmo) 

TFrDTRFCTi'^‘:),Ft,..23G 1 '^0 35 

DO 82 IIKslrTI 

82 IFtfiM.FU. lRt,Dilf I TK ) )G0 TO 35 , 

IFC tFLnwCftCO«ni‘' 3 F.O) . AnD . f STARl ( i''') . KO , FT MlSa ( iRKOiJf 1 1) 3 )•» 

1 GO TO 51 
35 COWTTNDE 

51 XF(FIjOW£KCOh,Nf'3-UT,TA-iT )lAMTsFOOw(Frrn%F;t,,) 

IFCIAMT.GE.O) GO TO ft 77 



60 


70 

o 

5 


TIsIT +1 

IRF!0f1(TI)alk;ij 

iF('FTNTSM('1i'.).-^t::,ST. -iff) GO TO <i5 
00 Y^'s T = 1,Tt 

TFI,ti'»*( TRPi;u?.i I r n = rf-’I;u'''M, THFl'tU T J ) 


Fiifiifl.' ( KOU'^ , Tt' ‘’iO’ r ! ) ) =:FLO 
It- .Of c-v I*"! 


Fon,. fisr j . 

ro*;TT-'‘YG'; 

■^t'T.O >U,K)=n 
RF.Ttj'^R 


CKCO' 

■■H 


• ta 


J R! I IM T ) 1-J 


' j sU 


2 B 





00100 
00200 
00300 
00400 
OObOO 
00000 
00700 
00800 
00900 
OlOOO 
01100 
01200 
01300 
01400 
01500 
01600 
01700 
OIBOO 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
02600 
02700 
02800 
02900 
03000 
03100 
03200 
03300 
03400 
1500 
1600 
1700 
1800 
03900 
04000 
04100 
04200 
04300 
044O0 
04500 
04600 
04700 
04800 
04900 
05000 
05100 
05200 
05300 
05400 
05500 
05600 
05700 
05800 


r 

c 

c 


c 

r* 


C 

c 

c 

c 


c 

c 

c 

c 

c 


c 


c 

c 


c 

c 

c 

c 




PROKRAMrifi* FO” THF r.^j^f^PAT jf'!,' 

OF f^FT-'flf-KS riO?.('='CTK.I 




0 I^'C 10(200 ,200) , R ' j-v:. ) 

T J’-eGt:R r,iA>Tt 150 ) ,Pt a TftP ( 3^0 } ,CA I-’ flR-) ) , T tf i' . 

I i.)0s2/ 

71^=21 

npF,.. (>'> TV = 1T0,F7L,P=:' , F rOr'-' , T.-i ' 1 

IPF: (M.* r ! = rou,»- TI.JTs'-Fxrh:*. .Ttifi'l 

1C Provides ttie sppo for rAoo'n,,. niri-her seneroi" iori 
done in the suoroijtinc 

n Ts'toe i.iMOer ot rando!^: nuooors nenerate'i oi 
one t.iT-e.if toe renuire of reo'’:iT. 
nu-r.oers Js Ji.ote tjien -t » toe anoroncf r.e 
FA'UjFH Ip aqair, exedtel. 

R Ts an array ot ra''S''io'n nutni-jers “feneratea. 


to toe netw 
in toe net.'' 
toe oronie.T 

tor any arc 


TC=36794 
Ns300 

CALL RANGEN(TC,.V,R,IFDFX') 

NNODE Is the number of no'les reaulrea 

NARC is the number of arccs required 

MCPMM is the number of commoflltles In 

MAXCAp Is the maximum canarltv allOAed 

DO 3001 Ni/MTlMssUb 
RSADCITN^f) N.vnDF,FAPC» f^’AXCAF 
DO 3001 MCnMMa?,5 

TvClD is the noc!e-node Incidence matrix 

DO 5 I = l,WNt)nE 
DQ 5 aal,N«!0DE 
TMCIDtr,J)=0 

Z'inE qlves the averaoe total dearee of arcs 
TOTDeg generates randomly the degree of ariv node 


Z>1ftF:=2*NARC/Mf<OUF 
00 35 T=l,MFnDF 
TOTDEGsRC I MDEX) =f2»7HPF 
INDEX=:IhDEX + l ^ 

IF tlFOFX.EO.iOCAIal. PAhGFN ( 1C , N , R # T wDEX J 
TFf iaTDEG.5T.2)itn:rEr:=2 


INDKC is the randomly generated Indegree of 
any node , ^ 

TOUOFG Is the randomly generated outdegree 
of any node 


TBlDEG=10Ti)FG»R C INUFX) 
lNDEXsI«DEX+l 



05900 

06000 

06100 

06200 

06300 

06400 

06500 

06600 

06700 

15 

06800 

06900 

07000 

07100 

c 

07200 

r 

07300 

r 

07400 

07500 

10 

07600 

07700 

07800 

07900 

12 

oaooo 

ORIOO 

082O0 

08300 

08400 

08500 

08600 

08700 

.25 

08800 

08900 

09000 

09100 

09200 

20 

09300 

09400 

09500 

40 

09600 

09700 

09800 

09900 

10000 

42 

10100 

10200 

10300 

10400 

10500 

lObOO 

10700 

15 

10800 

10900 

c 

11000 

c 

11100 

11200 

c 

11300 

c 

11400 

c 

H500 

c 

11600 

c 

11700 

c 

11800 

c 


I F .C r {sDEX . BO • 'OC^ 1.1. -4 f J C , ■ ' , t y ) 

IFCl^iKBG.LT.l ) l-.iCgrsJ 

fF (IhDBC.FO.liirOBGiTnrii'CsrT 4l)K^-1 

Tr)tlDBG5lOTnFG«T;.nBG 

,TFf.f .GB..'iCnMM)GO yn 75 

TP=0 

Dil 15 

TFd-ICTnfJ, U.CT.Ul Tf.rrt< + 1 

n:j 75 f-isl ,1 •■nfG“TP 

x r ■'•?.= u 


TCAX* Is the rhnho'nlv 'lehprate'i cap<?rit.li“s j^r 

^rcs store'l In rAn'io'^i.v nenprutet:^ 
locations of, l.he nooe-no'^e loctdenre ‘'^irjy 


Ks(>(Ti‘-nFX)*rr(V>i-: 

T'iOBy=Ti‘i!>KX + 1 

T F C If- OhX , BO . r. 3 C alt. P At'GFH (I.C , r! , P , T uhh' v ) 

TK(CK.Iit;.nCLP*!'n,f)H.{T.FO.Kl)flO TO 3o 

rcAP=H(i'Mjp:x)*PAyc:Ap 

TT^‘n,'.;=iTrtB+i 

r^'OKXsT'inEXfi 

XF (TiSiDFy.By.i'OCALT. PAdGF-J f I." , « , P , r'inEX ) 
IF(iCAP.BO.d)Gn T(1 17 

T F ( C I NC X n C K , ,r ) , wF, 0 ) . AFP, ( TT'IE , LF . FhOijF) 1 GO TO 30 

rftfCir)(K,T)=iCAP 

TMCID(I,K)=-ICAP 

CQNTIf^UE 

IF(I,EO.MiMOPE)GO TO 35 

rp=o 

no 20 J=l,MNODB 

IFCIMCIP(I,J) ,GT,01IP=Tp+i 

conunub 

no 35 W=l,IOUl)BG-lP 
ITME=0 

KsRC INDEX) ♦fJWODB 
INDEX=INDEX+1 

IF ( INDEX, FO, id CALL PANGFA; ( 1C , F # P , INDEX ) 

IF (CK.LE.FCnM’d.GP.d.EO.K)) GO TD 40 
ITMBsITMB+1 

icap=p(inuf:x)^faxcap 

TNnFXsINDEXtl 

IF ClWDEX.Ep.ldCALL PADGFN ( iC ,M ,P , InDEX ) 


IFClCAP.Fy^O) GO TO 42 


TFC(T'hCTPrK,T).tF,r)) 
JMCIf>(T,F) = lCAP 
TnCID(K,T)=-TCAP 

CiHTIl'jOfc 


AWI). CTTME.LP.Hi^nuE) )Gn TO 40 


from the node-node incidenre matrix ,,the input data 
Is qenerated 

K10DE Is the total number of nodes including the 
super source 

NUMARC Is the total numoer of. arcs includinq the 
arcs connecting the super source and source 
STARTCJ) Is thestartina node. of the Jth arc 
FINISH(J) Is the finishlno point of the Jth arc 
CAPCJ) is the capacity of the Jth arc 

3' 2 ' ' 



11900 

12000 

12100 

12200 

123^0 

12400 

12500 

126O0 

12700 

12800 

12900 

13000 

13100 

132O0 

13300 

13400 

13500 

13600 

13700 

13800 

13900 

14000 

14300 

14400 

14500 

14600 

14700 

14800 

14900 

15000 

15100 

15200 

15300 


15500 

15600 

15700 

15800 

15900 

16000 

16100 

16200 


!3J 111) I = i 

tto jal 

T ( T vC r D ( T , T K i,F . n ) Cii rn \ i 
. -<l;.mr3C = .V‘J.iAHC + 1 
31 -/TMi^Osr 

j.i,iA'7C) = r '.ic mm: r ,i!3 
110 cn'rit..'iF, 

00 I 2 r. 

3i:.’VR'!'{X) = lVf.i)^-,+ ) 

PI' ISrUDsT 
C6P( :r )=9')oo 

K OLP + t 

TFlr-'Or'Af .3r'(/)fin to iOOl , 

WRITS! Ton, *) i-nf-iAnr, 'U'nR'n 

WhTXS! T')1.1,T3 CSTAPT'(u) ,.7=1 ,7 U '-'APCI 
^IRITK! lun, t ) (FT,vjSH(,I) lArf ) 

wRiTSdon,*) fcftpCJ) ,J=i ,*'UMAr?c) 

3001 C1J''T10?,IE 

C005E (I(■^^1T=^I'I) 

CoOSF (IfMXTsTrjU) 

STC’P 

FID 

c <<<<<<<<<<<<<<<<<<<<<«<<«<<«<<<<<<<<<<<<<<«« 

c Tnls subrontlnp Is useti for tiir* 

c veneration of random noThers 

c >>>>>>>>>>>>>>>>>>>>>>>>>>>>>»>>>>»>>>>>>>>>>>> 

c Tfils routine uses the TPPf,. sunrom. ine HGOo 

c which Qonerates ranont.. numhers. 

SUBROUTINE RA WOfc.f' ! IC , , F , t NORX 1 

D1ME<<sTON RC3003 

IC=IC+387 

CAT.L GGUB(IC»^#P) 

I.'iDEXsJ 

RETURN 

EMD 



